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Abstract. Recycling strange stars to millisecond periods is studied within the framework of general relativity. We 
employ equations of state of strange quark matter based on the MIT Bag Model, with massive strange quarks 
and lowest order QCD interactions. The presence of the crust of normal matter is taken into account, with a 
bottom density assumed to be equal to the neutron-drip one. The calculations are performed by solving the exact 
2-D equations for rigidly rotating stationary configurations in general relativity. Evolutionary tracks of accreting 
strange stars are computed, and their dependence on the initial strange star mass and on the fraction of the 
angular momentum transferred to the star from the infalling matter, is studied. The differences between recycling 
strange stars and neutron stars are pointed out. 
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1. Introduction 

The problem of the spin up of the neutron (or strange) 
stars is an important element of the studies of recy- 
cled pulsars. It is widely believed that millisecond pul- 
sars {P < 10 ms) are old neutron stars (age 10^ y), 
with low magnetic field {B ^ 10^ G), which had been 
spun up by the accretion when they were members of the 
low-m ass binaries (see, e.g., van den Heuvel & Rappaport 
1992 ) . This evolutionary scenario of formation of millisec- 



ond pulsars recently obtained some observational support 
with the discovery of the first accretion-powered millisec- 
ond pulsar SAX J1808.4-3658 (Wijnands & van der KHs 
1998). 

The detailed theoretical modeling of the spin-up sce- 
nario requires the solution of several specific problems that 
involve the knowledge of the neutron star structure, that 
in turn depends on the equation of state of neutron star 
matter. In particular, one has to determine the accreted 
mass needed to spin up a compact star to the observed 
millisecond periods, within a reasonable time, in low-mass 
X-ray binaries. Finally, one would like to know what limits 
the maximum accreted mass - whether the upper bound 
results from the mass-shedding limit, or from some in- 
stability against gravitational collapse. The answers to all 
these questions depend on the (largely unknown) equation 
of state of dense matter at supranuclear density. 
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Kluzniak & Wagoner (1985) studied the spin up of 
an accreting neutron star (NS) for several equations of 
state of dense matter, using the lowest-order expansions in 
stellar angular momentum (slow rotation approximation, 
pioneered by Hartle 1967). Cook, Shapiro & Teukolsky 
( |1994b| , hereafter CST) modeled the recycling of NS with- 
out any restriction on neutron-star rotation, using their 
previous results for stationary rigidly rotating stellar con- 
figurations, which were obtained in the framework of gen- 
eral relativity for fourteen equations of state of dense mat- 
ter (Cook et al. 1994c ). A semianalytical model, ba sed on 
the CST results, has been used by Burderi et al. ( 1999 ) 
to study the dependence of the recycling scenario on the 
magnetic field. 

A hypothesis that the true ground state of matter 
is a deconfined self-bound quark plasma (Bodmer 1971, 
Witten 1984), called strange matter, implies the possibil- 
ity of the existence of strange quark stars (Witten 1984; 
first detailed models in Haensel et al. 1986, Alcock et al. 
1986). Recently, it has been argued that the accreting com- 
pact objects in some low-mass X-ray binaries could be 
strange stars (SS), and not NS (Bombaci 1997, Li et al. 
1999a,b). 

Previous studies of the properties of rotating strange 
stars pointed out significant, and in some cases strong and 
even qualitative differences between rapidly rotating SS 
and NS. In particular, the difference in the distribution 
of matter in the interior of the rapidly rotating NS and 
SS was shown to imply a significant difference in the outer 
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spacetime. Specifically, qualitative differences in the prop- 
erties of the innermost stable circular orbit (ISCO) around 
these two classes of rapidly rotating compact objects have 
been found (Stergioulas et al. 1999; Zdunik et al. 2000, 
Zdunik & Gourgoulhon 2001). As shown in Zdunik et al. 
(2001, hereafter Paper I), a careful treatment of the crust 
of SS to some extent softens these conclusions, but im- 
portant differences still remain. Recently, parameters of 
marginally stable orbits and properties of the accretion 
disk around SS have been studied, and compared with 
those around NS, by Bhattacharyya et al. (2001). 

Differences between rapid rotation of SS and NS, as 
well as those between properties of the space-time around 
them, indicate the need for a study of differences between 
recycling by accretion of SS and NS. In the present paper 
we study spin-up of SS by accretion of matter from the 
innermost stable circular orbit. We neglect the effect of 
the magnetic field, accretion, and radiation drag on the 
location of the ISCO, and on the infall of accreting mat- 
ter, which is justified for B ^ 10^ G and M <C MEdd- 
Our approach is thus similar to that of CST. We have re- 
cently calculated the models of rotating strange stars with 
crust in general relativity (Paper I) . Using the method de- 
scribed in Paper I, we determine the parameters necessary 
to calculate the motion of the matter accreting from the 
ISCO in the equatorial plane. The equation of state of 
strange matter, as well as the method for solving the rel- 
evant equations in general relativity, are briefly described 
in the first part of Sect. 2. The remaining part of Sect. 2 
is devoted to the formulation of the basic equations gov- 
erning the evolution of an accreting SS. Numerical results 
of computation of the evolutionary tracts of accreting SS 
are presented in Sect. 3, where we also study the basic 
differences between spin-up by accretion of SS and NS. 
Finally, Sect. 4 is devoted to the discussion of our results 
and gives some conclusions. 

2. Calculations 

The equation of state and calculational method for the 
stellar model are the same as that used in Paper I. Our 
EOS of strange matter, composed of massless u, d quarks, 
and massive s quarks, is based on the MIT Bag Model 
with the standard values of the Bag Model parameters 
for strange matter: bag constant B = 56 MeV/fm'^, mass 
of the strange quark nis — 200 MeV/c^, and QCD cou- 
pling constant ac — 0.2 (Farhi & Jaffe 1984, Haensel et 
al. 1986, Alcock et al. 1986). This EOS of strange quark 
matter (called SQMl in Paper I) had been also used in 
Zdunik et al. (2000) and Zdunik & Gourgoulhon ( ^001] ). It 
yields an energy per unit baryon number at zero pressure 
Eo = 918.8 MeV < E{^^Fe) = 930.4 MeV. The maximum 
allowable mass for static strange stars is M^^^. = 1.8 Mq. 
The crust in our model is described by the BPS EOS below 
the neutron drip (Baym et al. 1970). 

The general relativistic models of stationary rotat- 
ing strange stars have been calculated by means of the 
code developed by Gourgoulhon et al. (1999), which re- 



lies on the multi-domain spectral method introduced in 
Bonazzola et al. (1998). In the present calculation we use 
4 spatial domains (for details see Paper I). We supple- 
mented the version of the code used in Paper I with the 
equations determining the specific angular momentum Zms 
and specific energy Cms of a test particle falling from the 
ISCO to the surface of the star (see the Appendix). We 
determine also the angle /?, at which the particle falling 
freely from the ISCO hits the surface of the star, 



tan(/3) = 



EL 



(1) 



where p"", are radial and tangential momentum of the 
infalling test particle, in the frame comoving with the 
matter at the stellar surface. For energetic considerations 
(heating of the surface), as well as for determining the 
probability of crossing the electrostatic barrier at the sur- 
face of SS, the crucial parameter is also the energy of the 
particle as seen by the observer comoving with the star. 
This energy is defined by the time component of the par- 
ticle 4-momentum: E ~ . The equations defining the 
parameters of the particle hitting the surface of the star 
are presented in the Appendix. 

The evolutionary tracks of the accreting SS are calcu- 
lated using the simplest scenario in which the increase of 
the angular momentum of the star is related to the ac- 
creted rest mass increment by 

SJ SMb , (2) 

where we introduced the parameter xi - the fraction of 
the total angular momentum of the particle transferred to 
the star. In the idealized, limiting case the accreted mass 
element transfers all its angular momentum to the star 
{xi = 1). Equation (||) allows us to determine the total 
angular momentum J and and total rest (baryon) mass 
Mb, after accreting matter of rest mass AMb- 

Mb = Mb,o + AMb , (3) 

I-AMb 

J = ■h+ Xlln,siMB,J)dMB , (4) 

Jo 

where Mbo and Jq are initial baryon mass and angular mo- 
mentum of the star, respectively. The values of Mb and 
J determine then in a unique way the instantaneous con- 
figurations of accreting star, which form a one parameter 
family of rigidly rotating configurations, labeled by Mb- 
These configurations form evolutionary tracks of accreting 
SS. The evolutionary tracks will be labeled by the gravi- 
tational mass Mq of the initial configuration with baryon 
mass Mb,o- This initial configuration is assumed to be 
nonrotating, Jq = 0. 

For the star accreting at the rate Mb (as seen by a 
distant observer), the time evolution of J is given by the 
equation 

dJ , , > 

-7- = Xi Ims Mb , 

at 



(5) 



resulting directly from Eq. 
sured at infinity. 



where t is the time mea- 
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Fig. 1. Angular momentum and energy (per unit mass) 
of a particle at the marginally stable orbit along 
an evolutionary track for initial gravitational masses 
1, 1.4, 1.6, 1.72, 1.73, 1.75 M© and linear approximation in 
j = Jc/GM^ of Kluzniak and Wagoner (1985) (dashed). 
In the case of energy the curves for 1.73, 1.75 M0 cannot 
be graphically distinguished from the upper part of the 
1.72 Mq curve. The values of l^g and Cms for nonrotating 
configuration are ^/T2GM/c and 2\/2/3, respectively (ex- 
cept for the 1 M0 case, where the ISCO is defined by the 
stellar radius). 

3. Results 

3.1. Changes ofl^ns, Cms 3 long the evolutionary tracks 

As one can see from Eq.(H), the quantity determining the 
evolution of an accreting star is the specific angular mo- 
mentum of the particle at the ISCO: Ims- In the case of 
spherical, nonrotating configuration with R < Rms, the 
value obtained for the Schwarzschild metric is ^/l2GM/c. 
For a rotating star, inclusion of the lowest order correc- 
tion in a dimensionless parameter j — Jc/GM"^ lead s 
to the formula derived by Kluzniak & Wagoner ( 1985| ): 

\/T2(l — \ GM/c. Higher-order expansions in j can 

be found in Shibata & Sasaki ( |l998| ) (Eq. (B2) of their 
paper, where j is denoted by q). In this article, we use in- 
stead the exact formula for Z^s, as given in the Appendix. 

In Fig. |l| we present the value of the angular momen- 
tum and energy of the test particle at the ISCO along 
evolutionary tracks. The approximate formula underesti- 
mates the value of ^ms and the actual spin up of the star 
is larger than predicted by the linear approximation (in 
j) of Kluzniak & Wagoner (1985). The energy of the test 
particle at the ISCO, Cms, is smaller than that given by 
the linear approximation in j for high masses only (the 
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Fig. 2. The ratio /p"^, equal to the tangent of the in- 
cidence (impact) angle /3, at which a test particle falling 
freely from the marginally stable orbit hits the stellar sur- 
face. Thick solid line - exact results for nonrotating config- 
urations up to the maximum mass. Thin solid line - evolu- 
tionary tracks for Mq = I.6M0 and Mq = 1.7M0. Dashed 
lines - same quantities, calculated using the approximate 
formula (Kluzniak & Wagoner 1985). Evolutionary tracks 
terminate at the mass-shedding limit. 



initial mass Mg > 1.7 Mq). This is due to the fact that 
the next order terms in the expansion of Cms and /,„s in j 
are negative, while simultaneously rotational deformation 
of the star leads to the increase of both Cms and Zms- The 
negative correction, due to j, dominates only for the most 
massive stars, when the deformation due to the rotation 
is relatively small because of the very strong gravitational 
pull. Such a cancellation is not so pronounced for Cms, 
because there the relative magnitude of the p correction 



(see Shibata & Sasaki 1998 



is about twice that for l-^as 
Eqs. (B1,B2)). 

Figure |^ shows the value of the tangent of the inci- 
dence (hitting) angle of the accreted matter, falling freely 
from the marginally stable orbit, as measured in the rest 
frame comoving with the stellar surface. The deviat ion o f 
the approximate formula of Kluzniak and Wagoner (1985) 
from the exact results is quite large, reaching 50% at max- 
imum mass of the nonrotating configurations, and is even 
greater for rapidly rotating massive stars. The main rea- 
son for this discrepancy is the large gap between the stellar 
surface and marginally stable orbit for strange stars. For 
instance, at M = Mmax and i^iot = we have i? ~ 10 km 
and RrcLs — 16 km, so that the linear approximation of the 
expansion of tan/3 in (i?ms — R)/R, used by Kluzniak & 
Wagoner ( 198?: ), is unjustified. 

In Figure ^ we present the value of the energy of the 
inf ailing particle (minus its rest mass ttiqc^) as seen by the 
observer at the stellar surface. This energy can be inter- 
preted as the kinetic energy of the particle which should 
be lost as it gets a part of the star. As one can expect the 
energy tends to as we are approaching Keplerian limit 
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Fig. 3. Kinetic energy of a test particle falling freely from 
the marginally stable orbit as it hits the stellar surface 
(in the frame comoving with the surface of the star). 
The curves correspond to fixed rotational frequencies: 
v-sox. — 0, 1000, 1200 Hz. Thick lines - exact results (dashed 
part - ISCO defined by the stellar radius). Thin dotted 
lines - nonrelativistic approximation with relativistic time 
correction (see text). 

with R > Rjns (left part of v^ot = 1000 Hz curve), since 
then the Keplerian velocity at the equator is equal to the 
velocity of the rotation. The exact results are compared 
with the approximate Newtonian formula (with a correc- 
tion involving the relativistic time delay which changes 
the velocities as seen by the observer at the massive star) : 



SE = 



1 {Vrr.s-V? GM 

2 (1 - 2GM/Rc^) R 



1 



R 



R„ 



(6) 



where v = flR is the velocity of the star at the equa- 



tor and Vn 



^^ms^ms is the velocity of the particle at 



the marginally stable orbit. This formula approximates SE 
with a surprisingly small error. 

3.2. Accreted crust 

Standard models of strange matter predict the existence of 
an electron layer, which extends beyond the sharp quark 
surface by some ~ 300 — 400 fm. The electrons in this layer 
are bound to the quark surface by a positive Coulomb po- 
tential, produced by quarks, of height ~ 30 — 50 MeV. 
Spherical accretion of a plasma onto a SS with low mag- 
netic field {B <i 10^ G) was studied in (Haensel & Zdunik 
1991, Miralda-Escude et al. 1990). In the case of spheri- 
cally symmetric accretion, protons at the instant of hitting 
a bare (no crust) SS surface have an energy many times 
larger than the height of the Coulomb barrier, and tra- 
verse it with a large probability. After being absorbed by 
the strange matter, the protons (nucleons) dissolve into 
quarks in an exothermic process. 

As we see in Fig. H, in the case of accretion from the 
ISCO, nuclei in the infalling plasma hit the SS surface at 



Fig. 4. Gravitational mass versus central baryon num- 
ber density along several evolutionary tracks of accret- 
ing strange stars. Solid lines correspond to the typical 
case - all angular momentum of matter falling from the 
marginally stable orbit is transferred to the star (x; = 1). 
Dashed lines represent three cases for a Mq = 1.4 M© star 
{xi — 0.66, 0.55, 0.5 from the left to the right), and the case 
xi — 0.5 for a Afo — 1 M© star. Long-dashed, nearly ver- 
tical line at maximum mass is defined by the onset of in- 
stability with respect to axisymmetric perturbations (see 
the text), dotted line - Keplerian (mass shedding) limit. 



a very oblique angle f3 ^ 10°; they should therefore be ef- 
ficiently deflected by the electrostatic barrier. The kinetic 
energy of a proton hitting the SS surface can be estimated 
as — 300 MeV (Fig. ^) , and depends very sensitively on 
the value of fi. The largest value of 6E corresponds to 
the nonrotating star with maximum mass, since then the 
two components defining this energy are very large: the 
velocity of the freely falling particle relative to the nonro- 
tating star (the kinetic energy) and the gap between the 
marginally stable orbit and stellar surface (the potential 
energy). However a very small incidence angle increases 
strongly the effective thickness of the Coulomb barrier to 
be penetrated, as well as the effective thickness of the 
electron layer to be traversed (implying therefore higher 
energy loss). All in all, infalling protons (or more generally 
- nuclei) are not able to penetrate the Coulomb barrier, 
and one expects a gradual building-up of a normal-matter 
envelope (crust) on the top of the strange matter sur- 
face. The bottom density of the accreted crust is limited 
by the neutron-drip. Further accretion induces irreversible 
absorption of neutrons by strange matter, and their fu- 
sion with strange matter in a strong-interaction process 
n — > 2d + u, followed by the weak-interaction process 
u + d — > s + u. These reactions lead to an additional heat 
release (deconfinement heating) at the crust-strange mat- 
ter interface, of some 30 — 40 MeV per accreted nucleon. 
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3.3. SS parameters along evolutionary tracks 

In Fig. ^ we plot the gravitational mass of the star as a 
function of the central density, for different choices of Mb,o 
and for Jq — 0. This figure can be directly compared with 
Fig. 1 of CST, obtained for NS. The main differences be- 
tween neutron and strange stars spun up by the accretion 
are the following: 

— the difference between maximum masses for rotating 
and nonrotating SS is more than 0.4 Mq compared 
with at most 0.3 Mg for NS; 

— the specific angular momentum l^g for SS is sufficiently 
high to result in a characteristic "turning point" of the 
evolutionary track in the pc — M plane, at which the 
central density of accreting SS starts to decrease. The 
effect of the increasing angular momentum (leading to 
the deformation of the star and decrease of pc) starts 
there to be dominating over the the effect of the mass 
increase tending to increase the central density. 

— the rapid increase of J, compared to NS, determines 
the final fate of the accreting SS: nearly all evolution- 
ary sequences terminate at the mass-shedding limit. 
Instability with respect to axisymmetric perturbations 
(which would eventually lead to the collapse of accret- 
ing SS into a black hole) can be reached only for ini- 
tial (static) mass very close to the maximum one for 
xi = 1; for our EOS this is restricted to the narrow 
range 1.73 Mq < Mq < M^ax = 1-8 Mq. 

— in the case of accreting SS, there exists a region of 
stable rotating configurations with high M and high 
pc, which cannot be reached via accretion of matter 
onto an initially static SS (unless xi < 1). 

— contrary to the case of NS, for astrophysically inter- 
esting masses (say, M > Mq), we have, during accre- 
tion, i?oq < -Rms- This means, that along the evolution- 
ary track the marginally stable orbit is always located 
above the stellar surface. This conclusion is different 
from that presented in Paper I, because it refers here to 
the supramassive configurations, which have not been 
considered in our previous paper. 

The mass-radius relations for accreting SS along evo- 
lutionary tracks are shown in Fig. |^. The region that can- 
not be reached by the accretion is not so spectacular as in 
the case of Uc — M dependence. It is clearly visible that 
stars which, due to the accretion, directly enter the region 
of instability with respect to axisymmetric perturbations 
(implying collapse to black hole) lie in a very narrow re- 
gion (1.9 > M/Mq > 1.73, 11 km > i? > 10.2 km). 
However it should be noted that in principle it is possi- 
ble that stars spun up by the accretion to the rotational 
frequencies very close to the Keplerian one eventually are 
slowing down due to some processes leading to the angu- 
lar momentum loss (violating of course Eq. @j) and reach 
the region forbidden from the point of view of accretion. 

In Fig. P we present rotational frequency and grav- 
itational mass of the star - two quantities which could 
in principle be observed. It should be stressed that the 
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Fig. 5. Mass versus the (circumferential) equatorial ra- 
dius, along several evolutionary tracks of accreting strange 
stars. Notations as in Fig 0. 




Fig. 6. Frequency of the rotation versus gravitational 
mass along several evolutionary tracks of accreting strange 
stars. Notations as in Fig ^. 

point in the f — M plane do not uniquely defines the 
rotating configuration. In particular the crossing point 
between the Keplerian (dotted) and the instability limit 
(thick dashed) curves at about 2.5 Mq, 1400 Hz do not 
correspond to the same star (the Keplerian limit star has 
a total angular momentum J = 4.2 GMq/c and central 
density Uc ~ 0.5 fm~^ while the marginally stable star 
has J = 2.6 GMq/c and Hc = 1.1 fm~'^). As a conse- 
quence the region below the instability line and above the 
Keplerian limit (dotted) curve for supramassive stars can- 
not be treated as forbidden (stable rotating configurations 
exist there) although it can be reached by the accretion 
only in the case of xi < 1. 

The existence of the region which cannot be reached 
by the accretion can be suitably discussed in the J — Mb 
plane - two basic parameters from the point of view of 
accretion (Fig. |^). The evolutionary tracks are defined 
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Fig. 7. Angular momentum of the star versus its baryon 
mass, along several evolutionary tracks of accreting 
strange stars. At each point, the slope (first derivative) 
of evolutionary curves is equal to the value of xilms- 
Notations as in Fig. ^ 

by the requirement that the change of the stellar angu- 
lar momentum is connected with the accreted mass via 
Eq. (||), i.e., the derivative dJ/dMe along evolutionary 
tracks is equal to xilms- The allowed region is bounded by 
the Keplerian limit (high J) and by the instability with 
respect to axisymmetric perturbations (high Mb) lead- 
ing to the collapse. As it can be seen from Fig. ^ the ac- 
creting star eventually collapses if the evolutionary track 
crosses the instability limit (long dashed curve in Fig. 0), 
which is possible if d J/ dMe along this limit is larger than 
(dJ/dMB)evoi along the evolutionary track. For xi — 1 
this condition limits the possible parameters of collapsing 
accreting SS to a very narrow range (high Mb and low J) 
whereas for a smaller xi all configurations on the instabil- 
ity curve can be reached by accretion. It should be noted 
that, contrary to Figs. ^ and ||, Fig. ^ presents stellar 
configurations in a non-unique way, i.e., above the long- 
dashed curve there exist stable (lower ric, higher R) as well 
as unstable (higher ric, lower R) configurations. Stars be- 
low this curve do not exist, which is the direct consequence 
of this stability criterion (minimum J at fixed Mb). 

Figure ^ shows the rotational frequency of a spun-up 
SS versus the accreted baryon mass. The evolution of the 
rotational period depends critically on the fraction of the 
angular momentum at the ISCO, transferred to the star, 
Xi . On the contrary, the character of the spin- up depends 
rather weakly on the initial mass of the star. For xi = 1, 
the submillisecond periods are reached after accretion of 
about 0.3 M0. Accretion of 0.2 Mq is needed to spin up 
the SS to 1.56 ms (/ = 641 Hz), the minimum observed 
period of a millisecond pulsar. Notice, that at a fixed Mb, 
the abscissa in FigJ| is proportional to the accretion time, 

iacc = MB,acc/MB. 

An interesting quantity, in the astrophysical context, 
is the accreted rest mass required to spin up the SS to 
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Fig. 8. Rotational frequency of rotating strange star as 
a function of the accreted baryon mass. Dotted lines - 
Mo = 1,1.73,1.75 M0. Dashed lines correspond to xi = 
0.5 for Mo = 1 Mq and xi = 0.5 (dashed), xi = 0.55 
(dot-dashed), xi = 0.66 (long dash) for Mo = 1.4 M©. 
The last points of the each evolutionary track are labelled 
by the initial gravitational mass and denoted by a star 
(mass-shed limit) or circles (collapse). 

the termination point of the evolutionary track (i.e., to 
the mass-shedding limit, or to the threshold for instabil- 
ity with respect to the axisymmetric perturbations, which 
results in collapse into a black hole). For the evolution- 
ary tracks which eventually lead to collapse into a black 
hole, this maximum accreted mass depends sensitively on 
the initial mass Mo, i.e., how close Mo is to the maximum 
allowable mass. For example, for Mo — 1.75 M©, the max- 
imum accreted mass is only MB,acc = 0.07 Mq, while for 
the lowest Mq, for which the evolutionary track terminates 
by instability with respect to axially symmetric perturba- 
tions, Mq = 1.723 Mq, the SS accretes MB,acc = 0.2 Mq 
before collapsing into a black hole. The maximum accreted 
mass, required to spin up a SS from fiot = to the mass- 
shedding limit u-cot — i^K, practically does not depend on 
the initial mass (provided it is not too small, Mq > 1 Mq), 
and is roughly equal MB,acc = 0.6 Mq, i.e., is larger by 
some 50% than in the case of neutron stars (see CST). 

The parameters of the crust (the baryon mass and the 
thickness) are presented in Figs. |[ |l^. The interesting dif- 
ference with respect to the curves presented in Paper I for 
a fixed baryon mass of the star Mb is the slower increase 
(or even decrease) of Mb, crust and Ai?crust with increasing 
Vrot- This is caused by the increase of the SS mass dur- 
ing accretion: a larger mass corresponds to higher grav- 
itational pull and smaller crust (Paper I). Although the 
rotation results in the increase of Mb, crust and A_Rcrust at 
slow rotation rates the mass increase can dominate, lead- 
ing to the shrinking of the crust. Of course, the effect of 
the mass increase is more important, if the accreted mass 
is connected with a slower acceleration of the rotation (i.e. 
xi < 1). For example, for xi = 0.5 the size of the crust of 
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Fig. 9. Baryon mass of the crust along evolutionary 
tracks. Mq = 1, 1.4, 1.6, 1.72, 1.73, 1.75M0. Higher initial 
mass corresponds to lower Mb. crust- Dashed curves corre- 
spond to xi — 0.5 for Mq = IMq and xi = 0.5,0.55,0.66 
(from bottom to the top) for Mq — lA 




SS, and do not depend on the specific EOS of self-bound 
strange quark matter. 

The standard models of strange matter (a recent re- 
view can be found in Madsen 1999), predict the presence 
of a small fraction of electrons, which is due to the SU(3) 
symmetry breaking in the quark sector (m^ S> mu,md). 
The electrostatic potential, which bounds electrons to 
SS, produces then a huge outward-directed electric field, 
which is necessary to support a crust of normal matter 
on the strange matter core. Using such a strange matter 
model, one is able to study the formation and evolution 
of the crust during the accretion process. 

Very recently, Rajagopal & Wilczek (2001) demon- 
strated that quark matter in the color-flavor locked (CFL) 
superconducting phase of QCD is strictly neutral, despite 
unequal quark masses: n„ — nd — ng. If strange matter is 
actually in the CFL superconducting phase, then SS would 
be bare, with no electrons, and with a surface thickness 
of ~ he/ fi, where fj, ~ 400 MeV is the quark chemical po- 
tential. In absence of any repulsive Coulomb barrier, the 
stable coexistence of normal and quark matter would then 
be excluded. All accreted matter would be then absorbed 
in an exothermic process. However, the effect on the evo- 
lutionary tracks of accreting SS would be small. The gap 
Rins ~ Rcq would be larger, and the shedding-limit termi- 
nations would shift to somewhat higher frequencies, allow- 
ing for accreting of a larger amount of matter. However, 
the parameters of the marginally stable orbit, and the rate 
of spin-up of accreting SS would not change, because the 
maximum fraction of stellar mass and moment of inertia 
contained in a normal crust is at most 10^"*. 
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an accreting Mq — 1.4M0 SS star is continuously decreas- 
ing along the evolutionary track. 

4. Discussion and conclusion 

Our study of SS recycling by accretion have shown sig- 
nificant differences with respect to the case of NS. These 
differences stem both from the basic difference in the in- 
ternal structure (density profile) and from the difference in 
the outer space-time around rotating SS and NS. The lat- 
ter implies differences in the properties of the marginally 
stable orbit, from which the matter is falling onto the 
compact object, transferring its angular momentum to the 
star. Our calculations have been performed for a specific 
EOS of strange matter. However, the qualitative features 
of evolutionary tracks, described in Sect. 3, are generic for 



B^r^ sin' 
-r^d(p^) 



' 9{dip - N^dtf 



(A.l) 



where N, iV^, A and B are four metric functions of (r, 6). 
For a complete description of the coordinate system, met- 
ric functions and relativistic equations see Bonazzola et 
al. ( 1993| ) and Gourgoulhon et al. ( 1999| ). Here we present 
only the equations for the parameters of the particle falling 
freely from the marginally stable orbit (ISCO). 

The location (coordinate r^s) and orbital frequency 
rims of the ISCO are obtaine d by s olving the equations 
presented in e.g. Cook et al 1994a , Datta et al. 1998, 
Schaab & Weigel 1999 . For any particle falling along a 
geodesic from the ISCO to the stellar surface, the con- 
served quantities are the specific energy 



'Pt — Tms {N + Br,^,N'^Vn,s) 



(A.2) 
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and the specific angular momentum 



^ms ~ Pip r 



Br V 

3 ^ ' ms * ms 5 



(A.3) 



where denotes the covariant components with respect 
to the coordinates {t, r, 9, ip) of the particle 4-momentum 
(per unit mass), Kns is the particle 3- velocity as mea- 
sured by the zero-angular-momentum observer (ZAMO) 
and Fijis the corresponding Lorentz factor: 



B 



ms 



-1/2 



(A.4) 



(A.5) 



In the above equation, all the metric factors are to be 
evaluated at the location of the ISCO. 

Let us now compute the impact angle and energy of 
the particle as it reaches the stellar surface. First we con- 
struct an orthonormal tetrad connected with the observer 
comoving with surface of the star, the angular velocity of 
which is r2 (as seen from infinity). Let (ej-, e^, e^, e^) 
denote this tetrad. It is expressed in terms of the the co- 
ordinate basis (et, e^, er, eg) as follows: 



= u = — Bt + It — Bu, 

t J\f ^ AT ^ 



N 



F F 

^U^e, + — (1 
A* Br 

1 



A 
1 

Ar 



Br 



(A.6) 



eg 



where u is the 4-velocity of the stellar matter, Ug its 3- 
velocity with respect to the ZAMO, 



and F the corresponding Lorentz factor: 
F:= . 



(A.7) 



(A. 



It is easy to check that this tetrad is orthonormal, i.e. 

% ■ ^ ^ ^t^'^ ^ Diag(-1, 1, 1, 1). 

The freely falling particle that hits the surface of the 
star has a 4-momentum p with two conserved compo- 
nents: the quantities pt = — e,ns and = l^s given by 
Eqs. ( |A.2| ) and ( [A.3| ). For a motion within the equatorial 
plane, the component pg vanishes. The component Pr is 
obtained from the normalization relation p • p = — 1: 



pI = a^ 



V^ms * ''ins } 'n 



iV2 



52? 



- 1 



(A.9) 



which of course is nonzero at the surface of the star unless 
the ISCO is determined by the stellar radius. 



In the frame comoving with matter (tetrad A.6) the 
momentum of the particle is given by: 



P 



(en 



p • 6' 

Br 



1 



A2 

pA%Y 

BrUs 



-pti%y 



(A.IO) 



N 



The angle /3 at which the particle falling freely from 
the marginally stable orbit hits the surface of the star is 
given by 



P 

tanp — — 

p-p 



(A.ll) 



The energy release per unit mass due of infalling mat- 
ter in the comoving frame is p* — 1, i.e. 



N/T] 



(A.12) 



This is the value which should be used for the energy 
balance considerations at the stellar surface (heating of 
the matter). If we assume that all of this energy is emitted, 
this results in a luminosity (per unit accreted mass) seen 
by a distant observer SE/u* = e,ns ~ ^Ims ~ N/T, as it 
has recently been demonstrated by Sibgatulin & Sunyaev 

( pooo| ). 
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